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DYNAMICAL LOCALIZATION OF DIRAC PARTICLES IN 
ELECTROMAGNETIC FIELDS WITH DOMINATING 
MAGNETIC POTENTIALS 

JEAN-MARIE BARBAROUX, JOSEF MEHRINGER, EDGARDO STOGKMEYER, 

AND AMAL TAARABT 


Abstract. We consider two-dimensional massless Dirac operators in a radi¬ 
ally symmetric electromagnetic field. In this case the fields may be described 
by one-dimensional electric and magnetic potentials V and A. We show dy¬ 
namical localization in the regime when lim IRI/IAI < 1, where dense point 

r—>-oo 

spectrum occurs. 


1. Introduction 


Graphene is a two dimensional material consisting of carbon atoms arranged in a 
honeycomb lattice which was isolated in 2004 [14] . Behind its remarkable properties 
such as Klein tunneling and finite minimal conductivity |9] stays the fact that at 
low excitations energies the dynamics of charge carriers is described by the massless 
two-dimensional Dirac operator |3|. For technological devices based on graphene 
one needs the ability to confine and control the mobility of charge carriers. However, 
confining Dirac particles is not an easy task due the so-called Klein effect, where 
particles are able to penetrate electric potential walls |9] with very little reflexion 
index. In |7] it was argued that in presence of rotational symmetric electric and 
magnetic fields one could confine or deconfine Dirac particles by manipulating the 
strength of the fields at infinity, i.e., far away from the sample. Our main result is a 
dynamical statement on this effect and a continuation of a recent work m by two 
of the present authors. Before presenting the result let us explain this phenomenon 
with more mathematical details. 

Denote by H the two-dimensional massless Dirac operator coupled to a radially 
symmetric field E = Ef on the plane and a radially symmetric transversal magnetic 
field B. If the fields are sufficiently regular H is a self-adjoint operator densely 
defined in L^(]R^,C^) acting as 


( 1 ) 


H = cr- {-iV - A) -f K, 


where V : —>■ M is the electric potential satisfying 


K(x) 



E{s)ds = V (r) 
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(abusing notation, we write V{r) to denote b^(x), where r = |x| £ [0,oo) is the 
standard radial variable). Here cr = {ai,a 2 ) is a matrix-valued vector whose com¬ 
ponents are the first two Pauli matrices 



The vector potential A = (^ 1 ,^ 2 ) : —>■ generates the magnetic field B, with 

B = ^ 1^2 — 92^1- We choose the rotational gauge, i.e., we set 

A(x) := -A{r) ( , with A{r) = - f B{s)sds. 

We note that, besides some local regularity requirements for (V, A), H is essentially 
self-adjoint on independently of the growth rate of A and V (see [4]). 

In this setting there exists a unitary transform [18] (see also [101 Section 6]) 

( 2 ) —> L^(R+,C^)^e^{Z), 

such that the operator H can be written as a direct sum of operators on the half-line 

(3) 

iez 


where 

(4) hj = —ia 2 dr + (Ti{A—^) + V on L^(R''',C^), 

with rrij = j + ^ for j £ Z. Clearly, the spectra of H and hj are related through 

(5) a{H) = (J a{hj), adH) = [J and <7pp{H) = [J (Jpp{hj). 

iGZ jez 


For the operators hj we have the following properties, assuming sufficiently regular 
fields: If A{r) —>• 00 as r —>■ 00 and 


( 6 ) 


lim 


Vir) 


A{r) 


< 1, 


then the spectrum of hj is discrete for each j £ Z (see [121 Proposition 1] for the 
precise regularity conditions). In contrast, if H(r) —>■ 00 as r —>■ 00 and 


(7) 


lim 

r—¥oo 


Air) 


Vir) 


< 1, 


then the spectrum of hj equals the whole real line and it is purely absolutely 
continuous [HI Propsition 2]. This suggests delocalized particles in the regime given 
by 0 and confined particles in the one given by 0. However, the latter is not 
obvious since, for fields satisfying 0, H may have dense point spectrum (see [13] 
and m Theorem 7.10] for the case when B decays at infinity and [12] for the case 
when B is not decaying at infinity). We recall that dense point spectrum may lead 
to non-trivial dynamics. In fact, in this case, it is only known that the wave-packet 
spreading is sub-ballistic [17] (the result stated in m is for Laplace-type operators 
but can easily be adapted for the Dirac case). Moreover, there are examples of 
systems with pure point spectrum where the spreading rate is arbitrarily close to 
the ballistic one 016] (see also 0)- 
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Concerning dynamical results we know that particles in electromagnetic fields 
satisfying © behave ballistically, i.e., for any finite energy state ip G 
and K > 0 one has (see M) 




POO 


Jo 



dt - for large T > 0. 


The main result of this work is that under condition ([6]) the operator H exhibits 
dynamical localization, i.e, for any /c > 0, for any finite energy interval I, and 
for any state ip G T^(R^,C^), with sufficient regularity in the angular variable 
(depending on k; see (fT^ l. holds 


sup 

t>o 




2 

< OO, 


where Pi{H) is the spectral projection of H onto I. 

Let us now state our assumptions and result more precisely: 


Hypothesis 1. A^V & 


(8) 

\A{r) 

—>■ OO 

(9) 

lim 

V{r) 


r—>00 

A[r) 

(10) 

lim 

A'ir) 


r—>oo 

A^{r) 


and they satisfy 

as r —>■ OO, 

< 1, 

= 0 . 


Recall that U (see @) is the unitary map that decomposes H in the direct sum 
of the operators hj. For a given ip G L'^ (R.^, ) we write 

(11) Uip = 0 Lpj, with G L^(R~'’,C^). 


Our main result is the following theorem. 


Theorem 1.1. Let k > 0, / C R &e o bounded energy interval and let Pi{H) be 
the spectral projection of PI onto I. Assume that A and V fulfill Hypothesis{J\ and 
let Ip G L^(R^,C^) be a normalized state such thatlAip = © satisfies 


(12) 

IjT WTjf < OO. 


jez 

Then we have 


(13) 

sup X Pi{H)ip < OO 


t>o 


We note that the condition HU) is related to regularity of the initial state ip 
in the angular variable. Indeed, let r~^l'^ip with ip G L^(R''' x [0,27r),C^) be 
equal to ip G L^(R^,C^) expressed in polar coordinates. Then, m follows if 
iP G iJ'"/2([0, 27r), L^{R+)) © iJ«/2([0, 27r), L^{R+)). Here, for k > 0 

([0, 27r), L2 (r+)) := g 27t) x R+), ^(1 + K|)« \\ue\\l2^R+) < oo}, 

iei, 


is the fractional Sobolev space on the torus [8] and ui is the £ — th Fourier coefficient 
of u with respect to the variable 9. To be more explicit, note that Pip = © pj where, 
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for g e L2([o,27r),C2), 

{Tg){ 3 ) := ^ Me and Mg = 


.ie/2 


We notice that 


(-^5)0') = 


\i-92ij + 1 )/ 

The angular momentum operator J := —ide+cfzl2 satisfies TJT* 
for simplicity that n = 2n with n S N, we have 

)'/" = ||j"^|| = ||Ha, + ^r^ 


0 ^ 

i^-^e/2 J ■ 


rrij. Assuming 


< 



{-ide)^ 


< T' 


V' 


ff"([0,27r),L2(K+))2 ■ 


We can also derive a sufficient condition on for m to hold. To avoid unnecessary 
complications, we stick to the case of even values of k. Since de = —X 2 dxi +xidx 2 , 
we have 


{de'N 


2 


{-X2dj;^ + xidx2)'''fp\\'^ < Ck ^ ||(1 + |xp)2i9“i/’ 

|a|<fc 


where for the multi-index a = ( 01 , 02 ) we used the notation 9“ = d', 
in Cartesian coordinates, condition (IT^ holds if ^ belongs to the weighted Sobolev 


Hence, 


space 


,C2) := {gGL^R^,^^) \ (1-f |x|2)f g(x) e C^)}. 


2. Proof of the main result 


The strategy of the proof is as follows. We first rewrite the moment in ((T51) using 
the decomposition of i? as a direct sum of hj ’s and the representation 1 /; ~ 0 
(see (O). The main idea consists in taking advantage of the discrete spectrum of 
the operator hj and the exponential decay of its eigenfunctions to compensate the 
growth of the moment. However, the decay is not uniform in j, it turns out (see 
(1^ in Lemma HTTl) that the exponential decay estimate can only be derived outside 
an interval which grows with j. This is due to the fact that the term A — mj/r 
in hj can be controlled by a fraction of A only outside a ball of radius Xj with 
rj\A{rj)\ ^ |mj|. Hence, we split the space K'*' x Z into two regions through the 
function /, given in Definition 12.21 In the first region, where the values of r|A(r)| 
are sufficiently large compared to those of j, we control the growth of the moment 
in r with the exponential decay of the eigenfunctions obtained in Section |T] by an 
Agmon-type argument. We also need to control the number of these eigenfunctions 
in the spectral region we consider. This latter bound is derived in Section [3] by 
using arguments due to Bargmann. In the second region where the values of j are 
large compared to r|A(r)|, the regularity assumption (1121) in the angular variable 
for our initial state yields a decay in the variable j which is used to control the 
growth of 

We now turn to the detailed proof of Theorem 11.11 

Since Uip = we have 


( 14 ) 
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We note that it is enough to consider only the sum for |j| > Jqj for sufficiently 
large Jq > 1. Indeed, let Nj be the number of eigenvalues of hj in the interval 
(k) 

I. Let , k = 1,2,Nj denote the corresponding eigenfunctions of hj. Then, 
expanding Pi{hj)(pj in terms of the we have 


E 

bi<Jo 




< 


Nj 

EE| 

\j\<Jok=l 




2 

< OO. 


where the last inequality holds in view of Remark 14.21 below. To estimate the right 

_ nOO 

hand side of (HI we split the integration / into two regions characterized by 

3 

(15) Tj = rj{5o) := sup{r G R+ | \mj \ > dor|A(r)|}, 

for some Sq G (0,1) which is chosen in the proof of Lemma [4.II below. 


Remark 2.1. Note that since A is continuous and A{r) oo as r ^ oo we have 
that Tj < OO and that Vj —>■ oo as |j| —>■ oo. Moreover, the supremum in (TTSl) is 
attained and hence 


\mj \ = 6orj\A{rj)\. 


\mj\ < dor\A{r)\, r>rj. 


(16) 

Moreover, we note that 
(17) 

Definition 2.2. Let 6 G C'°°(K.''', [0,1]) with 9{r) — 0 for r < 1 and Oir) = 1 for 
r>2. We set f{r,j) = fj{r) := 0{r/Mj) and f^ := 1 - fj. 

Proof of Theorem \l.li We have 

2 


E 


L-«/2g -ithj 


Pi{hj)ipj 


(18) 


< 2 ^( f,r^l\-^^^^Pi{h,)^, + 

jGZ ^ 

We first estimate the second term of the right hand side of (fT5|) using the regularity 
in the angular variable for the initial state as given by In what follows we 

pick Jo > 1 so large that |A(r)| > 1 for all r > rj and |j| > Jq. In particular, we 
have that 

(19) rj < \mj\/5o. 

Using this and the support properties of /°, we have 


( 20 ) 


E |/ 

bi>'^o 


^ E 

\f^{6r,r^^e-^*>^^Pj{h,)^, ' 


bi>-^o 


< ^ (6|mj|/(5o)''||e Pi{hj)\\^ \\ipj\\^ < oo , 

\3\>Jo 

where we used m in the last bound. 
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We now estimate the first term in the right hand side of ((T51) . For |_;/| > Jq we 
compute 

2 

/ W -• \ 

2 


< sup ^ 
11011=1 \ ^ 


, fe=i 


Ni 


Ni 


- XI ^ sup ^ 

l.--, " 11011 = 1 


(0 ^\|2 


N, 


1^1 


=i: 


/2f(fc) 


fc=l 


Consider the function fj defined at the beginning of Section S] below. Since fj = 
fjfj we have, choosing also Jq > J2 (see Lemma IFTI) 




< 








(fc) 




< 


«;/2p-7e 




Ti 


where p{r) = |A(s)|ds is the exponential weight defined in Lemma o Since 
r'^/'^Q-'^e decays monotonically at infinity, we may choose Jq > 1 to be so large that 
the supremum of is bounded above by Hence 

< .^(3^^.)'‘/2g-7(e(3rj)-e(2rj))^ 




Note that due to (HB 


^ 37-7 


g{3rj) - g{2rj) = I r|H(r)|— > -^Ind) > 1^, 

J2rj I’ CO 0Oo 

for |j| large enough. Thus, using that Vj < \mj \ and Lemma lOI we get for |j| > J 2 

k^l ^ 

2\mj I 

< S'^C^C/Imj-r+Mnlmjle”'^^. 

Since the last bound is summable for |j| = — 1/2| > Jq we get the expected 


result. 


□ 


3. Estimate on the number of eigenvalues of hj 

Let r be a self-adjoint operator on a Hilbert space H with purely discrete spec¬ 
trum. We set for an interval / C K 

Ni{T) := dim Pj{T)n, 

i.e. Nj{T) denotes the number of eigenvalues of T in / counted with multiplicity. 

Lemma 3.1 (Bound on the number of eigenvalues for hj). There is a Ji > 1 such 
that for any E > 0 there is a constant Ce > 0 so that 

(21) N[_E,E]{hj) < CE|m_;-|ln|mj|, for \j\ > Ji. 
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Proof. We first note that 

hj = -i(T 2 dr + tTi (A(r) - ■PV{r) 

is essentially self-adjoint on (7“(R+,C^). In addition, we obtain by the spectral 
theorem 

In the sense of quadratic forms on we obtain for any e G (0,1) the 

estimate 

>(I - e)[-ia2dr - cn^ + (JiA{r)Y + (l - \)V‘^{r) 

=(1 - e)[( - +A‘^{r) - ^V^(r) - asA'ir) - '^A{r)\. 

Let i5 := (1 — e)/2. Due to (l9|) and (fTOjl we have that 5A^ — a^A' and eA? — I/^/e 
are positive outside a large ball i? for e G (0,1) sufficiently close to 1. Let C := 
\\V‘^/e + a^A'Wi^oatgy Then we find 


h]>{l-e) ( - i(J2dr - 0-1^)^ -b SA^{r) - C - ^^\A{r)\ 


We write 


where 


/i| - > (1 - e) ( - ia 2 dr - -b Wj 


Wjir) = SA^ir) - C - (1 - - ^^\A{r)\. 


Let Rj := rj{6/4), where rj{-) is defined in (ITSl) . This yields, |TOjj < |r|yl(r)| for all 
r > Rj. Moreover, we may pick Ji so large that (recall that Rj —>■ oo as |jj —>■ oo) 

^A‘^{r)-C-{l-e)-^E'^ >0, for all r>Rj,\j\>Ji. 

Thus , Wjli^R.^^) > 0 and 

( 22 ) - L;' > (1 - e) [( - ia2dr - ■ 

Define 

(23) Dr.= {rGiO,Rj) \ |m,| > |r|bl(r)|}. 

Note that if r G (0, Rj) fl (K'*' \ Dj) then |bl(r)^ — |^(r)| > 0. Hence we have 

(24) hi - £;2 > (1 - e) [( - ia2dr - a+ W<], 
where 

(25) W<{t) := (M2(r) - M|H(r)|)l,,, - (C^ + (1 - e)-'L;")l(o,ij,]. 

An application of the min-max principle leads to 

N[o,e^] [h]) = -/V(-oo.o] [h] - E'^) 

< y-ooM{{ - -<Jl^y+W<). 

A direct computation shows that 

( - ia2dr - CTi^)^ = -d1 + ^mArrij - erg) 

-d1 + ^ rrijiruj - 1) 0 

0 -dl + :^mj{mj + l)^ 
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Note that mj{mj ± 1) > 0 for |j| > Ji. Using the generalized Bargmann estimate 
[2] (see also Theorem XIIL9]) we obtain for \mj\ > 1/2 

-^(- 00 , 0 ] ( “ ^ mj{mj — 1 ) + 

^ 2^^ ir ^1^/ (^) if ^3 > h 

^ Wld?- ifTOj<-i 


< 


. 2 |m ,-1 + 1 Jo 


and 


-^(- 00 , 0 ] { ~ dr ^ rujinij + 1 ) + +/') 

'' 2^+1 ^ 

. fo°" r\W<{r)\dr if <-\ 


< 


and therefore 


(26) iV(-oo.O](( - i <^ 2 dr - CT1+)' + W <) < 

Now we estimate using the definition of Dj 


1 


- 1/2 


7 ’|+^(r)|dr. 


(27) 


> ^ (C'+(l-e)-lU 2 )i ?2 /• , , 2 lrr,-l 

r|lU^(r)|dr < --- -+ / r\5A^{r) -Lr^|A(r)||dr 

2 J Di 


< 


(C'+(l-e)-iU 2 )i ?2 


2 |TOjP(r)|dr. 


'Di 


Furthermore, 


/ 2|TOj||A(r)|dr = / 2|mj||A(r)|dr + / 2|mj ||A(r)|dr 

(28) dDj Jo^inio,!) JDjnii,oo} 

<2|m,|||+^,.[o_i] + ^ln(i?,). 

Note that in view of Remark l2.ll and the fact that |^(?’)| grows at infinity, we have 
for sufficiently large Ji 

4|mj| ^ 4 |toj| 


mR,)\ 


111 > Ji- 


This together with (1251) . (|271) . and (|25)) yields the result. 

4. Exponential decay of eigenfunctions of hi 


□ 


Let rj = rj{So) be given as in (fTK)) . We note that do G (0,1) will be fixed through¬ 
out the proof of the next lemma. For the function 6 as defined in Definition 12.21 
we set f{r,j) = fj{r) := 0{rlr^). 

Lemma 4.1. There exist 7 > 0 and J 2 > 1 such that for all |j| > J 2 the following 
holds: Let G L^(R+,C^) be a normalized eigenfunction of hj with energy E G I. 
Then, for some C > 0 (independent of j), 


(29) 




~ r. 


where for r > 0, g{r) := / |A(s)|ds. 

Jo 
















Remark 4.2. It is clear from the proof that the exponential decay of the eigen¬ 
functions of hj remains true for |j| < J 2 , however, in this case we get a different 
constant in front of the exponential. 

Remark 4.3. Throughout the proof of Lemma \4-.1\ we use that hj and kj := hj — 
V are essentially self-adjoint operators on (R"'", ) (see for instance [a and 

references therein). Moreover, we also use that V is a perturbation with respect to 
the magnetic Dirac operator kj in the sense that D{V) D and there 

exists C such that 

IlFv^ll < C{\\kM\ + IIV^II) for all e C'o°°(R+, C^). 

Indeed, let ip G C'“(R“'",C^) and XR be a smooth characteristic function of a ball 
of radius R> 0. We set Xr = 1 — XR o-nd V = where := Vxr and 

:= Vxr- We thus have 

||V"¥>|| < < Cr II^II + Ux^rpW , 

since is bounded and V is dominated by A at infinity according to assumption 
Moreover, for R large enough, we use cni) and the identity k'j = —df -\- Aj — 
(73 A' to write 

5 WAx'kTW < WkjXRTW < Wi^x'kM + WXRkjTW < llVXfill ||‘/?ll + \\kjp\\ ■ 

Proof of Lemma \4.1\ In order to derive the Agmon-type estimates, we follow [lOj . 
We set 

(30) A, :=(A-^), 

and we notice that \Aj \ > (1 — ^o)|A| on the support of fj. Let fj be a normalized 
eigenfunction of hj associated to an energy E. We define 

(31) g, :=e^^'/^-, 

Q 

where 7 S (0,1) and = - - such that gir) = / |A(s)|ds. Note that g^ is 

1 + eg Jo 

bounded and differentiable. Consider the operator 

(32) kj = —ia2dr + aiAj = hj — V, 
and we define 

(33) Qj := Re{kje'^‘^‘gj, kjQ-''^^^g,). 

The task is to obtain bounds for Qj that will allow us to bound g uniformly in e. 
Lower bound. Notice that 

(34) [kj,e^^‘] = -ia2XQ'eA<^% 
so that we rewrite 

Qj = Re(e-'^^' g^, g^) 

= ^e{{kj - ixg)a2)gj, (kj + i 7 g> 2 )fo) = \\kjgjf - 7 ^ • 

Moreover, we have 

(35) k'j = —d(: + Aj — 0-3 A'. 

In view of Remark [Q and m for any e > 0 there exists Jg > 0 such that for all 
|j| > Jg one has 
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and therefore 


(36) 


to. - ^3A'j)gj) > ((1 - Sof - e){gj,A^gj). 


Now we drop the term —d^ of (1351) . This together with (p6)) yields 
Qj>{{l-Sor-~e)\\Ag,f-j^\\g',g,f 

> ((1 - do)^ - e) \\Ag,f - 7 ^ ||pto|l' = ((1 - - e' - 7 ^) \\Ag, 

Upper bound. We rewrite 

Qj = Re{kje'^‘^‘gj,fj{E - V)ij) + Re{kje^‘^‘g^, -icr2fj^j) 

= Re{e^^^g,J,iE - + Re(e^^'g„ [k,J,{E - U)] Q 

+ Re{e<^^gj,-ia 2 fj{E - V)^^) + Re(e^‘?'gj, kj, -*0-2/) ^j) 

= Re{e'*^‘gjJj{E - + Re(e^^'gj, [%, -icra/j] ^j), 


Since 


Re(e^^'5,, -*a2/'(U - V)^,) = Re(e^^'5,, ^ kj] fj^j) = 0. 

Furthermore, we use 

Re(e^'^'ff„ /,(U - ¥)%) = ||(U - U) g,f . 


In addition, we find some C > 0 such that 

kj,-i<y2fj] ^i)| = i-fj + 2cr3to/j)G)l 




< C- 


p7e(2r'i 


^\\9j\\ + \\A,g,\\) 


”T) 

—^(lto-|| + (i + toP5.ll), 
u 

where in the last inequality we use that rj > 1 (for sufficiently large |j|) and (l30l) 
together with the support properties of fj. We thus get 

p-fgi2rj) 

(38) Q, < IKU - U) 5,f + C-dig,II + (1 + do) |tog,||). 

U 

Then, combining (|37)) and (p 8 l) we get for |j| > Jg 

p7e(2r,) 

(39) {g^,{{{l-Sor-i-^^)A^-{E-Vf)g,)<C -(|| 5 ^.||+ (l + do) Pg,||). 

U 

According to ([H]) and Remark 12.11 we may pick dp, e and 7 so small that there are 
constants Jso,e,'y,csg,e,'Y > 0 such that |A| > 1 on the support of fj and, for all 

IjI > 

(40) (g„ [((1 - do)^ - e~- 7^)A2 - (E - V^jg,) > csg,,.^ Ug.f ■ 


This together with (IMl) yields 


(41) ||Ag,||< 




(ltoll/Pff.ll + (i + to)< 


C eieC^Ti) 


(2 + do). 


C<5o,e,7 C(5o,e,7 j 

The claim follows using the theorem of monotonic convergence for the limit e —>■ 0 
of dH]). □ 
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